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Abstract: 

Different proposals for the wave function of the universe are analyzed, with an 
emphasis on various forms of the tunneling proposal. The issues discussed include 
the equivalence of the Lorentzian path integral and outgoing - wave proposals, the 
definitions of the outgoing waves and of superspace boundaries, topology change and 
the corresponding modification of the Wheeler - De Witt equation. Also discussed 
are the "generic" boundary condition and the third quantization approach. 
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1. Introduction 



In quantum cosmology the whole universe is treated quantum-mechanically and is de- 
scribed by a wave function rather than by a classical spacetime. This quantum approach to 
cosmology was originated by DeWitt more than 25 years ago, and after a somewhat slow 
start has attracted much interest during the last decade. The picture that has emerged from 
this line of development 2^,4,5,6,7,8,9,10 j g a sma rj dosed universe can spontaneously nu- 
cleate out of nothing, where "nothing" refers to the absence of not only matter, but also of 
space and time. 

The wave function of the universe is defined on superspace, which is the space of all 
3- metrics hij(x) and matter field configurations <f>(x), 



It is invariant under 3-dimensional diffeomorphisms and satisfies the Wheeler-DeWitt equa- 
tion [1] 



Here, Ti is a second-order differential operator in superspace. In principle, ip(h, 4>) should 
contain the answers to all meaningful, questions one can ask about the universe. However, 
the conditions necessary to specify the appropriate solution of eq.(1.2) and the procedure by 
which information can be extracted from that solution are far from being understood. 

As (almost) any differential equation, the Wheeler-DeWitt equation has an infinite num- 
ber of solutions. To get a unique solution, one has to specify some boundary conditions in 
superspace. In ordinary quantum mechanics, the boundary conditions for the wave function 
are determined by the physical setup external to the system under consideration. In quan- 
tum cosmology, there is nothing external to the universe, and it appears that a boundary 
condition should be added to eq. (1.2) as an independent physical law. 






(1.2) 
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Several candidates for this law of boundary condition have been proposed. Hartle and 
Hawking 7 suggested that ip(h, <p) should be given by a Euclidean path integral over compact 
4-geometries g^(x,r) bounded by the 3-geometry h^j{x) with the field configuration <f>(x), 

r{h,4>) 

ij = J [^][#]exp[-^(^»]. (1.3) 

In this path-integral representation, the boundary condition corresponds to specifying the 
class of histories integrated over in eq.(1.3). Compact 4-geometries can be thought of as 
histories interpolating between a point ("nothing") and a finite 3-geometry hij. 

A Euclidean rotation of the time axis, t — > —ir, is often used in quantum field theory 
because it improves the convergence of the path integrals. However, in quantum gravity the 
situation is the opposite. The gravitational part of the Euclidean action Spj is unbounded 
from below, and the integral (1.3) is badly divergent. Attempts to fix this problem by analytic 
continuation 11 were only partly successful, and at present it remains unclear whether one 
can meaningfully define an integral such as (1.3). 

Alternatively, I proposed 10 ' 12 that ip(h, <p) should be obtained by integrating over Lorentzian 
histories interpolating between a vanishing 3-geometry and (h, 0) and lying to the past of 

(M), 

lKM)= / [dg][d<t>]e lS . (1.4) 
This wave function is closely related to Teitelboim's causal propagator 13 ' 14 K(h,2, 02|^1> </>i), 

ip(h,<p) = K(h,m- (i-5) 

Linde 8 suggested that, instead of the standard Euclidean rotation t — > — ir, the action 
Se in (1.3) should be obtained by rotating in the opposite sense, t — > +ir. This gives a 
convergent path integral for the scale factor, which is all one needs in the simplest min- 
isuperspace models. But in models including matter degrees of freedom or inhomogeneous 
modes of the metric one gets a divergent integral. Additional contour rotations might fix 
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this problem, but no specific proposals have yet been formulated. Halliwell and Hartle 15 
discussed a path integral over complex metrics which are not necessarily purely Lorentzian 
or purely Euclidean. This encompases all of the above proposals and opens new possibilities. 
However, the space of complex metrics is very large, and no obvious choice of the integration 
contour suggests itself as the preferred one. 

In addition to these path-integral no-boundary proposals, one candidate law of boundary 
conditions has been formulated directly as a boundary condition in superspace. This is the so- 
called tunneling boundary condition 16 ' 17 which requires that ip should include only outgoing 
waves at boundaries of superspace. The main weakness of this proposal is that "outgoing 
waves" and the "boundary of superspace" have not been rigorously defined. The Lorentzian 
path-integral proposal (1.4) was originally suggested 10 as a path integral version of the 
tunneling boundary condition, and indeed the two proposals give the same wave function 
in the simplest minisuperspace model 18 . In the general case, the equivalence of the two 
proposals is far from being obvious. 

I should also mention a completely different approach to quantum cosmology, the so- 
called third quantization 19 ' 20 ' 21 ' 22 ' 23 ' 24 . Here, the wave function of the universe ip is pro- 
moted to a quantum field operator and is expressed in terms of creation and annihilation 
operators for the universes. The problem of defining the boundary conditions is then re- 
placed by the problem of determining the in-state of the quantum field vp. With the radius 
of the universe playing the role of time, it is argued that creation of universes from nothing 
corresponds to an "in- vacuum" state at vanishing radius. 

The problem of boundary conditions for the cosmological wave function is related to 
the problem of topology change in quantum gravity. In the path integral approach, one has 
to specify whether the integration in (1.3), (1.4) is performed over 4- manifolds of arbitrary 
topology, or only a restricted class of topologies is included. In the tunneling approach, part 
of superspace boundary corresponds to boundaries between different topological sectors, and 
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one has to decide what kind of boundary condition should be imposed there. Moreover, we 
shall see in Section 7 that topology change not only affects the boundary conditions for ip, 
but also leads to a modification of the Wheeler-DeWitt equation. 

In this paper I shall review the status of the tunneling wave function of the universe and 
attempt a more precise formulation of the tunneling boundary condition. As a prototype 
for this boundary condition, the next section discusses the process of bubble nucleation in a 
false vacuum, which is in many ways analogous to the nucleation of universes. The outgoing- 
wave boundary condition for a nucleating bubble will be formulated using a spherical min- 
isuperspace model. In Section 3, similar approach is applied to the simplest cosmo logical 
minisuperspace model: a Robert son- Walker universe with a cosmological constant, A > 0. 
In Section 4, the wave function for the same model is obtained by analytic continuation from 
the "bound-state" wave function for A < 0. Section 5 discusses the Lorentzian path integral 
approach and its equivalence to the outgoing wave boundary condition. A possible extension 
of these approaches beyond minisuperspace is discussed in Section 6. There, it is suggested 
that some general properties of the potential term in the Wheeler-DeWitt equation may 
allow one to define outgoing waves in the general case. 

The issues of topology change is tackled in Section 7. It is argued that topology-changing 
transitions can occur through superspace boundaries, but generally involve configurations 
in the interior of superspace. This implies that the Wheeler-DeWitt equation needs to be 
modified. A possible form of this modified equation is suggested. Section 8 gives some critical 
comments on the third-quantization approach to topology change, and Section 9 contains 
some concluding remarks. 

2. Bubble Nucleation 

To discuss the nucleation of true vacuum bubbles in a metastable false vacuum, we shall 
make a number of simplifying assumptions. First, we shall assume that the bubble radius 
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at nucleation is much greater than the thickness of the bubble wall, so that the bubble 
can be approximated as an infinitely thin sheet. Second, we shall use the semiclassical 
approximation, assuming that the tunneling action is large. (This is always true for a thin- 
wall bubble, provided the theory is weakly interacting). The nucleating bubble is then nearly 
spherical and can be adequately described by a minisuperspace model with a single degree 
of freedom, the bubble radius R. Finally, we shall disregard the gravitational effects of the 
false vacuum and assume the spacetime to be Minkowskian. 

In our minisuperspace model, the worldsheet of the bubble wall is described by a single 
function R(t), and the Lagrangian is 

L = -4iraR 2 (l-R 2 ) 1 / 2 + ^-eR 3 . (2.1) 

Here, a is the wall tension, e is the difference between the energy densities of the false and 
true vacuum, and R = dR/dt. The momentum conjugate to the variable R is 

p R = 47raR 2 R(l - R 2 )- 1 / 2 (2.2) 

and the Hamiltonian is 

H=[p 2 R + (4naR 2 ) 2 ] 1 / 2 -^eR 3 . (2.3) 

Bubble nucleation does not change the energy of the system, and if the false vacuum 
energy is set equal to zero, we have 

H = 0, (2.4) 

which can be rewritten using (2.3), 

p 2 R + U(R) = 0, (2.5) 

U(R) = (A7raR 2 ) 2 (l - R 2 /Rq) , (2.6) 

where Rq = 3<r/e. The equation of motion for R(t) can be obtained from (2.2), (2.5) and 
(2.6), 

R 2 R 2 = R 2 -Rl, (2.7) 



and the solution is 

R(t) = (B$ + t 2 ) 1/2 - (2-8) 
The worldsheet metric of the bubble is 

ds 2 = (1 -R 2 )dt 2 -R 2 {t)dQ 2 , (2.9) 
where dfl 2 is the metric on a unit sphere. With a new time coordinate, 

r = R sinh -1 {t/R ) , (2.10) 
we recognize it as the metric of a (2+l)-dimensional de Sitter space, 

ds 2 = dr 2 - R 2 {T)dtt 2 , 

R(t) = R cosh(r/i? ) • (2.11) 

If the bubble wall gets inhabited by some 2-dimensional creatures, they will find themselves 
living in an expanding inflationary universe. If they are smart enough, they may also figure 
out that their universe was spontaneously created at r = 0, and thus eq.(2.11) applies only 
for r > 0. 

How would these 2-dimensional physicists describe the quantum nucleation of the uni- 
verse? In quantum theory, the energy conservation (2.4) gets replaced by 

W> = 0, (2.12) 

where ip(R) is the "wave function of the universe" and the momentum operator is pr = 
—id/dR. The square root in (2.3) is complicated to deal with, and it is much easier to use 
the energy conservation law in the form (2.5), 

[-d 2 R + U(R)]ip = . (2.13) 

The transition from (2.12) to (2.13) involves commutation of the non-commuting operators 
R and pr, which is justified, as long as 

\RpRiI>\»\[R,P R W\ = \ii>\, (2-14) 
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that is, away from the classical turning points, where pr ~ 0. Using the classical equations 
of motion for R(t), we find that (2.14) is violated in a small neighborhood of the turning 
point Rq, 

5R/Rq ~ (aRlr 2 « 1 . (2.15) 

Since the correct operator ordering is not known, we shall keep the simplest choice as in 
(2.13). 

We now come to the problem of determining the boundary conditions for ip{R). Only one 
non-trivial condition is required; the second would simply determine the overall multiplicative 
constant. The WKB solutions of eq.(2.13) for R> Rq are 

i/>±(R) =p(R)- 1 / 2 exp i^zi J p(R')dR' T^/^j , (2.16) 

where 

p(R) = [-U(R)} 1 / 2 (2.17) 
is the classical momentum. To the leading order in the WKB approximation, 

p R 4>±(R) « ±p(R)4>±(R) , (2.18) 

wherep^ = —id/dR. This shows that tp+(R) and ip-(R) describe, respectively, the expand- 
ing and contracting bubbles. In the quantum nucleation process, only an expanding bubble 
must be present, and thus we require that for R> Rq the wave function should include only 
the outgoing wave, if)+(R). 

In the classically forbidden range, < R < Rq, the two solutions of (2.13) are 

$± = b(i?)r 1/2 exp (±J*° \p(R)\dRj . (2.19) 

With the outgoing wave boundary condition at large R, the wave function in this range is 
determined 25 by matching at R Rq, 

il>(R < Rq) = ^+{R) + l -4>-(R) ■ (2-20) 
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The two terms on the right-hand side of (2.20) have comparable magnitude at R Rq, but 



A different choice of operator ordering would not affect (2.21), but it could affect the pre- 
exponential coefficient. 

Having obtained the result (2.21), the 2-dimensional physicist could be puzzled about 
its meaning What does it mean to find the nucleation probability for a bubble when there is 
only one bubble? Even if we assume that there are other bubbles, they are unobservable, so 
how can we test this theory observationally? Of course, in the case of a nucleating bubble, 
there is an external observer for whom the nucleation probability has a well-defined meaning. 
This may or may not be so in the case of the universe. But the point I want to make is 
that even a worldsheet observer can derive some useful information from the wave function 
of the universe. If, for example, several different types of bubble can nucleate, with different 
values of a and e, then the observer is more likely to find herself in the type of bubble with 
the highest nucleation probability (assuming, of course, that such bubbles are suitable for 



Furthermore, nucleating bubbles are not exactly spherical, and one could in principle 
calculate the amplitude for a bubble to have a given shape. This problem has been solved 
in the perturbative superspace approximation which includes all the degrees of freedom of 
the bubble, but treats all but radial motions as small perturbations 29 ' 30 . It turns out that 
perturbations of a spherical bubble can be represented as excitations of a scalar field $ that 
lives on the bubble worldsheet and has a tachyonic mass, m = -3i2§. The mode expansion 
of this field contains four "zero modes" which represent overall space and time translations 
of the bubble, while the remaining modes describe deviations from spherical shape. As in 



in the most of the forbidden range the ip+(R) term dominates. The exponential factor in 
the tunneling probability can be determined 2 6> 2 7,28 from 




(2.21) 



2-d life). 
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the cosmological case" 3 ' , one finds that the bubble nucleates with the field $ in a de 
Sitter-invariant quantum state. 33 This prediction should be testable both by external and 
worldsheet observers. 

Extension of this analysis beyond perturbative superspace is a very complicated problem 
which has not yet been solved. The bubble worldsheet can, in general, be represented in 
a parametric form as x^(£ a ) with a = 0,1,2. An external observer would evaluate the 
amplitude to find a bubble in a given configuration at x° = T by evaluating the path 
integral 



Given that there was no bubble at x = 0, the integration should be taken over all compact 
worldsheets bounded by the given 2-surface at x {) = T and satisfying < x Q (C,) < T. As I 
said, calculating the integral (2.22), or even making it well defined, is a very difficult problem. 

Fot a worldsheet observer, £° is a time coordinate and is a set of four interacting 

scalar fields. She would find the restriction on the range of unnatural and would 

probably define the no-boundary wave function ip in (2.22) as an unrestricted integral over 
x^(£). The two wave functions will generally be different, but in the semiclassical regime 
the integral (2.22) is dominated by the neighborhood of the classical path, and the wave 
functions will be essentially the same. It would be interesting to further investigate this 
connection between the bubble wave functions from worldsheet and target space points of 
view. At present, eq.(2.22) is purely formal, and its connection to the standard Euclidean 
formalism 27 ' 28 for calculating the vacuum decay rate is obscure. 

On a qualitative level, one expects quantum fluctuations to grow large at small length 
scales, and if large deformations are allowed, then the bubble wall can cross itself, and small 
"daughter bubbles" can be chopped off. When viewed at very small scales, the bubble wall 
may in fact have a fractal structure, with a dense foam of small bubbles surrounding it. 
Moreover, the worldsheet observer may discover that on sufficiently small scales her bubble 




(2.22) 
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is not a 2-d surface after all, but is more adequately described by certain solutions of (3+1)- 
dimensional field equations. Similar problems may face human observers as they explore 
distances approaching the Planck scale. 



3. de Sitter Minisuperspace 

Turning now to the cosmological wave function, we first consider the simplest minisu- 
perspace model, 

where p v is a constant vacuum energy and the universe is assumed to be homogeneous, 
isotropic, and closed: 

ds 2 = a 2 [N 2 {t)dt 2 -a 2 {t)<Kll}. (3.2) 

Here, N(t) is an arbitrary lapse function, dVL^ is the metric on a unit 3-sphere, and a 2 = 
2G/37T is a normalizing factor chosen for later convenience. Substituting (3.2) into (3.1), we 
obtain the Lagrangian 

r / 2 \ 

C = ±N[a(l-jp)- Aa 3 J , (3.3) 

and the momentum 

Pa = -ad/N, (3.4) 
where A = (4G/3) 2 p v . The Lagrangian (3.3) can also be expressed in the canonical form, 

C=p a d-NH, (3.5) 

where 

H = ~ (^ + a-Aa 3 ^ . (3.6) 

Variation with respect to p a recovers eq.(3.4), and variation with respect to N gives the 
constraint 

H = 0. (3.7) 
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The corresponding equation of motion for a is (for N — 1) 



a 2 + 1 - Aa 2 = , 



(3.8) 



and its solution is the de Sitter space, 



a(t) = iy _1 cosh(ift) , 



(3.9) 



where H = A 1 / 2 . 

Quantization of this model amounts to replacing p a — > —id/ da and imposing the Wheeler- 
De Witt equation 



and the parameter 7 represents the ambiguity in the ordering of non-commuting operators a 
and p a - This equation is very similar to eq.(2.13) for a nucleating bubble, and the following 
discussion closely parallels that in Sec. 2. 

The 7-dependent term in (3.10) does not affect the wave function in the semiclassical 
regime. Without this term, the equation has the form of a one-dimensional Schrodinger 
equation for a "particle" described by a coordinate a(t), having zero energy and moving in 
a potential U(a). The classically allowed region is a > if -1 , and the WKB solutions of 
eq.(3.10) in this region are 



7 + lTa- U(a) *' (0) =°- 



(3.10) 



where 



U(a) =a 2 (l- Aa 2 



(3.11) 




(3.12) 



where p(a) = [-U(a)] 1 / 2 . The under-barrier a < H 1 solutions are 




(3.13) 



For a » if -1 , 



Pa^±(a) ~ ±p(a)il>±(a) , 



(3.14) 
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and eq.(3.4) tells us that ip-(a) and ip=(a) describe an expanding and a contracting universe, 
respectively (assuming that N > 0). 

In the tunneling picture, it is assumed that the universe originated at a small size and 
then expanded to its present, large size. This means that the component of the wave function 
describing a universe contracting from infinitely large size should be absent: 

ip(a > H~ l ) =^-(a). (3.15) 

The under-barrier wave function is found from the WKB connection formula, 

^{a<H- l )=^ + (a)- l -^(a). (3.16) 

Away from the classical turning point a = if -1 , the first term in (3.16) dominates, and the 
nucleation probability can be approximated as^' 1 ^ 

-H- 



exp ( -2 jT \p{a')\da'\ = exp ("g^") ■ (3-17) 



It should be noted that the choice of N > in (3.4) is a matter of convenience. With 
the opposite choice, the roles of ip+(a) and ip-(a) would be reversed, and the boundary 
condition (3.15) would be replaced by ip{a > H^ 1 ) = i[)+(a). This would result in a time 
reversal transformation ip(a) ip*(a). Another way to look at this is to note that the time 
coordinate t is an arbitrary label in general relativity, and it is a matter of convention to 
choose time growing or decreasing towards the future (where "future" is defined, e.g., by 
the growth of entropy or by the expansion of the universe). Clearly, there is no physical 
ambiguity here, and once the convention is set, the tunneling wave function in this model is 
uniquely defined. 

At this point, I would like to mention the "generic" boundary condition suggested by 
Strominger 34 . He argued that since the nucleation of the universe is governed by small-scale 
physics, the boundary condition on ip should be imposed at small a, rather than at large a 
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as in the tunneling approach. The large-scale behavior of ip can then be determined without 
specifying the precise form of this boundary condition. The under-barrier wave function is 
generally given by a linear combination of ip+(a) and ip-(a), and for a "generic" boundary 
condition at a = 0, one expects the two terms to be comparable at small a. However, ip+(a) 
decreases exponentially with a, while V>— (a) exponentially grows and therefore dominates 
for all but very small a. The corresponding wave function in the classically allowed range is 
found with the aid of the WKB connection formula: 

iP(a < H- 1 ) =^_(a), 

i)(a>H- 1 ) = ijj + (a) + ijj-(a) . (3.18) 

The same wave function is obtained 35 by applying the Hartle-Hawking prescription to this 
model. 36 

The case for imposing boundary conditions at small a appears to me unconvincing. The 
same argument could be applied to bubble nucleation, but there we know that the correct 
boundary condition is the outgoing wave at large radii. Another familiar case when physics 
is confined to small scales while the boundary conditions are imposed at infinity is a bound 
state, like the hydrogen atom. In the next section, we shall discuss how the tunneling wave 
function can be obtained by analytic continuation from a "bound-state" universe. 

4. Tunneling Wave Function by Analytic Continuation 

The quantum-mechanical wave function for the decay of a metastable state is often 
obtained by analytically continuing the bound state wave function from the parameter values 
for which the corresponding state is stable. A similar approach can be adopted in quantum 
cosmology. As an example, we again consider the minisuperspace model (3.1), but now with 
p v < 0. In this case A < 0, and it is clear that the classical equation of motion (3.8) has 
no solutions. However, microscopic, Planck-size universes could still pop out and collapse as 
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quantum fluctuations. Then one expects the wave function to be peaked at very small scales 
and to vanish at a — > oo. 

When dealing with analytic continuation, approximate solutions like (3.12), (3.13) are 
not sufficient, since the neglected terms can become large after continuation. We shall, 
therefore, use the exact solutions to eq.(3.10) which can be obtained 37 for a particular 
choice of the factor-ordering parameter, 7 = — 1. With the boundary condition 

^(a^oo) = 0, (4.1) 

the solution is the Airy function 

tf(a) = Ai(z) , (4.2) 

where 

z = (-2A)- 2 / 3 (l - Aa 2 ) . (4.3) 
The asymptotic behavior of (4.2) at large a is 

1 / 2 exp[-(-A) 1 /2 a 3 /3 ^ (44) 

Continuation to positive values of A amounts to changing (— 2A)~ 2 / 3 — > (2A)~ 2 / 3 exp(=)=27ri/3), 
where the sign depends on the direction of rotation in the complex A-plane. Choosing the 
upper sign and using the relation 3 ^ 

2e ±«/3 Ai(ze T2 7 r l /3 ) = M ( z j ± ^ (4>5) 

we conclude that the wave function for A > is 

ip(a) = Ai(z) + iBi(z) , (4.6) 

with 

z = (2A)- 2 / 3 (l-Aa 2 ). (4.7) 
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This is the tunneling wave function 17 . The corresponding asymptotic form at large a is 

if>(a) oc a' 1 / 2 exp(-iA 1 / 2 a 3 /3) . (4.8) 

At this point, I would like to comment on one important difference between the above 
analysis and the standard treatment of the decay of a metastable state. In the standard 
approach, the Schrodinger equation for the bound state of a particle, 

Hip = Eif;, (4.9) 

is solved with the boundary conditions xjj — > at both x — > oo and x — > — oo. The energy 
eigenvalues E n are then completely determined by the Hamiltonian H = —d 2 + U(x). In 
the course of analytic continuation, as the parameters of the potential U(x) are changed, 
the eigenvalues E n also change and develop imaginary parts as the corresponding states 
become metastable. The resulting wave functions describe a probability that is exponentially 
decreasing with time inside the potential well by gradually leaking to infinity. On the other 
hand, in the quantum-cosmological model (3.10) the eigenvalue of the Wheeler-De Witt 
operator is fixed at E=0. At the same time, the wave function is defined on a half-line 
a > 0, and the boundary condition (4.1) is imposed only at a — > oo. The wave function is 
time-independent, and a steady probability flux at a — > oo is sustained by an incoming flux 
through the boundary at a = 0. In fact, as eq.(1.5) suggests, the tunneling wave function is 
more appropriately thought of as a Green's function with a source at a = 0, rather than an 
eigenstate of the Wheeler-De Witt operator. This will be further discussed in Sec. 5, 7. 

We note finally that a "generic" choice of boundary condition at a = would lead, for 
A < 0, to a wave function which is not confined to small scales, but instead increases without 
bound at a — > oo. 



5. Tunneling Wave Function from a Path Integral 
To discuss the relation between the outgoing-wave and path-integral forms of the tun- 
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neling proposal, we shall consider a slightly more complicated minisuperspace model: a 
Robertson- Walker universe with a homogeneous scalar field. After appropriate rescalings of 
the scalar field (f> and scale factor a, the corresponding Lagrangian and Hamiltonian can be 
written as 

C = ±ie a + e 5a (-a 2 + <p 2 -V(<j ) )], (5.1) 
n = l[e- 3a (-pl+pl) ~e a + e 3c V(0)] . (5.2) 

Here, a — In a, V((f>) is the scalar field potential, and the lapse function has been set 
N = 1. 

The path integral (1.4) for this model can be expressed in the form 14 ' 18 

roc 

K(q 2 ,qi)= / dTk(q 2 ,qy,T), (5.3) 
JO 

k(Q2,qi;T)=J [dq]exp^ij^ Cdt^j , (5.4) 

where q = (a,(f>) and the integration is taken over all paths a(t),<f>(t) beginning at q\ = 
(«!, (pi) at t — and ending at q 2 = (a 2 , (fi 2 ) at t = T. The function k(q 2 , q\,T) in eq. (5.4) 
has the familiar form of an amplitude for a "particle" to propagate from q\ to q 2 in time T 
and satisfies the Schrodinger equation 

^-n 2 ^k(q 2 , qi ;T) = (5.5) 

with the initial condition 

k(q2,qr,0) = 5(q 2 ,qi) . (5.6) 
The equation for K(q 2) q\) follows from (5.5), (5.6): 

H 2 K(q 2 , qi) = -iS(q 2 , q\) . (5.7) 

Here, H is the Wheeler-DeWitt operator, 

n = \e- ia \dl- d\-U{aM (5.8) 
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with "superpotential" 

U(a,<f)) =e 4a [l -e 2a V(<j))}, (5.9) 

and I am ignoring the factor-ordering ambiguity. The subscript "2" of H in (5.5) and (5.7) 
indicates that a and in (5.8) are taken to be a2 an d 02- 

Apart from an overall factor, the operator H in (5.8) is just the Klein-Gordon operator 
for a relativistic "particle" in a (l+l)-dimensional "spacetime", with playing the role of 
a spatial coordinate and a the role of time. The "particle" moves in an external potential 
U(a,(p). Let us now consider the behavior of K(q2,q\) as «2 — > ±00 with a\ fixed. We 
must first note that for a — > —00 the potentital (5.9) vanishes, and K(q2,q\) should be 
given 39 by a superposition of plane waves, exp[ik(a2 ±02)]- Since the path integral in (5.4) 
is taken over paths originating at some finite {a\, and going off to large negative 0:2 , this 
superposition should include only waves with k > 0. (Recall that p a > corresponds to 
a < 0). 

As «2 — *■ +00, the potential U(a,<f>) diverges, and the WKB approximation becomes 
increasingly accurate. The dependence of K(q2,q\) on q2 is then given by a superposition 
of terms e 1 ^ , where S is a solution of the Hamilton- Jacobi equation, 

2 _ 2 

da J \ d<p / 



In each term, the function S(a, <f>) describes a congruence of classical paths with 

d<f) (dS/d<p) 



(5.11) 



da (dS/dct) ' 

For V((f)) > 0, U(a,4>) w -e 6a \/(0) < 0, and it follows from (5.11) that |d^/da| < 1. 
Hence, the "particle" trajectories are asymptotically "timelike" and correspond either to 
expanding universes with p a = OS/ da < or to universes contracting from an infinite size 
with OS /da > 0. Since all paths originate at a\ < 00, the superposition should include 
only terms with dS/da < 0. For V{4>) < 0, the trajectories are asymptotically "spacelike" 
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and cannot extend to timelike infinity i+ or to null infinity I+. One expects, therefore, that 
K{q.2-> Ql) ~^ for q2 at i+ or 2+. 

Thus we see that the propagator K(q2,qi) satisfies the outgoing-wave boundary condi- 
tions both at a — > — oo and a — > +00. The tunneling wave function (1.4) is obtained by 
letting a\ — > —00 and integrating over all initial values of <j>, 

/oo 
d(f>'K(a, 0| - 00, 0') . (5.12) 
-00 

The trajectories then originate at the past timelike infinity i—, but the behavior of ip on 
the rest of the superspace boundary should be the same as that of K. This is illustrated 
in Fig. 1 for the case of V{4>) > 0. The probability flux is injected into superspace at i- 
and exits in the form of outgoing waves through J_ and i+. We conclude that the path- 
integral and the out going- wave forms of the tunneling wave function are equivalent, at least 
in the simple model (5.1). This is not very surprising, since eqs.(5.3)-(5.7) coincide with 
the standard equations for Feynman propagator in the proper-time representation, and the 
causal boundary conditions for the propagator are the same as the outgoing-wave boundary 
conditions for ip. 

6. Beyond Minisuperspace 

The main difficulty in formulating the outgoing-wave boundary condition in the general 
case is similar to the difficulty with the definition of positive-frequency modes in a general 
curved spacetime. There is, however, a hopeful sign. Our definition of outgoing waves in the 
minisuperspace model (5.1) was based on rather general properties of the potential U(a, 0): 
its unbounded growth at a — > +00 and its vanishing at a — > —00. It is not difficult to verify 
that the superpotential in the Wheeler-DeWitt equation has similar properties in the general 
case. 
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The general form of the Wheeler-DeWitt equation can be written as 1 

(V 2 -W = 0, (6.1) 



where 



V 2 = / d 5 xN 



ijkl 5hij8h k i 2 5<p 2 



(6.2) 



is the superspace Laplacian, 



X 1 icy 

G ijkl = 7, h ~ ( h ik h jl + h il h jk - h ijhkl) (6-3) 

is the superspace metric, 

U= I d 3 xNh 1 / 2 [-R^ + ]-h i U,i<t>,j + V((l>)} (6.4) 



is the superpotential, hij(x) and N(x) are, respectively, the 3-metric and the lapse function 
in the (3+1) decomposition of spacetime, 

ds 2 = (N 2 + N.N^dt 2 - 2N t dx l dt - h ij dx i dx j , (6.5) 

h = det(hij) and is the curvature of 3-space. As before, matter fields are represented 
by a single scalar field and I have ignored the factor-ordering problem. The wave function 
ip is a function of hij(x) and <f>(x), but is independent of N(x). The metric can be 
represented as 

= e 2a hij , (6.6) 

where det(hij) = 1. Then the Laplacian term in (6.1) is oc exp(— 3a), the first two terms 
in the superpotential (6.4) are oc exp(a), and the last term is oc exp(3a). The relative 
magnitude of these terms for a — > ±oo is the same as in eq. (5.9). 

In the limit a — > -co, the superpotential U in (6.1) becomes negligible, and one can 
hope to define outgoing modes analogous to the plane waves of the previous section. This 
possibility has also been suggested by Wald 40 in a different context. Here, I will not attempt 
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to analyze the most general case and illustrate the idea in a reduced superspace model which 
includes all degrees of freedom of the scalar field 0, but only one gravitational variable a. 
With the scalar field represented as 

<j ) (x) = (2n 2 ) 1 / 2 J2fnQn(x), (6.7) 

n 

where Q n (x) are the harmonics on a 3-sphere, the superspace Laplacian (6.2) takes the 
form 31 ' 32 




The plane-wave asymptotic solutions are then 

^(«, fn) = exp(ik a a + i ^ k nfn) (6.9) 

n 

with 

k 2 a -J2 k n = °- ( 6 - 10 ) 

n 

The tunneling wave function includes only terms with k a > 0. This is the boundary condition 
at a — > — oo. 

To formulate the tunneling condition on the remainder of superspace boundary, one first 
has to specify what that boundary is. In other words, we should decide what class of metrics 
and matter fields should be included in superspace. The form of the Wheeler-DeWitt equa- 
tion (6.1)-(6.4) suggests that we should include all configurations {hij(x), <f>(x)} for which 
^1/2^(3) ^ ^l/2^ij0 ^ an( j ^1/2^(0) are integrable functions. Then the superpotential U 
is finite everywhere in superspace and will generically diverge towards the boundary. This 
happens, in particular, at a — > +oo. As \U\ — > oo, some components of the gradient VS in 
the Hamilton- Jacobi equation 

(VS) 2 + U = (6.11) 

should also diverge, and one can hope that the WKB approximation will become asymp- 
totically exact, thus allowing one to define outgoing waves. 41 ' 42 For example, when some 
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dimensions of the universe become very large (e.g., a — > oo), the classical description of the 
corresponding degrees of freedom becomes increasingly accurate. Denoting these classical 
variables by q; and the remaining variables by qj, the asymptotic form of the wave function 
can be written as 31,43,44,45 

^(c,q) = J2 elSN{c) XN(c,q). (6.12) 
N 

The Hamilton- Jacobi functions S]^(c) describe congruences of classical paths, = —dS/dc^. 
The tunneling boundary condition selects the solutions of (6.10) which include only outgoing 
paths, evolving towards the boundary. 

It should be noted that superspace defined by the condition \U\ < oo includes a very 
wide class of configurations. The metric and matter fields have to be continuous, but not nec- 
essarily differentiable. In particular, scalar fields with discontinuous derivatives and metrics 
with 5-function curvature singularities on surfaces, lines, and points are acceptable configu- 
rations. This conclusion fits well with the path integral approach, where it is known that the 
path integral is dominated by the paths which are continuous, but not differentiable. The 
superspace configurations can be thought of as slices of these paths. 

If the definition of outgoing waves along the lines indicated in this section is indeed 
possible, then the same argument as in Section 5 suggests that the wave function defined by 
the path integral (1.4) should satisfy the outgoing- wave boundary condition. One advantage 
of the path-integral definition is that it may be consistent even if outgoing waves cannot 
always be defined. Another advantage is that the path integral version appears to be better 
suited to handle topology change (see Sec. 7). 

7. Topology change 

In the discussion, so far, I have not touched upon the issue of topology change in quantum 
gravity. This issue, however, can hardly be avoided, since the "creation of a universe from 
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nothing" is an example of topology-changing event. 

The Wheeler-DeWitt equation (6.1) is based on canonical quantum gravity, which as- 
sumes the spacetime to be a manifold of topology i?xE, where R is a real line and E is a 
closed 3-manifold of arbitrary but fixed topology. The corresponding superspace Qj^ includes 
only 3- metrics of topology E. We can define the extended superspace Q including all possible 
topologies. It can be split into topological sectors, with all metrics in each sector having the 
same topology. 

The division of superspace into topological sectors can be illustrated by lower-dimensional 
examples. In the (l+l)-dimensional case, 3-geometries are replaced by lines (strings), and 
topological sectors can simply be labeled by the occupation number of closed strings. In 
(2+1) dimensions, a point g G Q corresponds to a number of closed surfaces (membranes), 
and each surface can be characterized by the number of handles. 46 Each topological sector 
of Q can thus be labeled by an infinite set of integers {riQ,ni, ...} giving respectively the 
occupation numbers for surfaces with 0, 1, ... handles. In the (3+l)-dimensional case, there 
is a much richer structure, but a topological classification of 3-dimensional manifolds has not 
yet been given. 

Creation of a universe from nothing described in Sections 3-5 is a transition from the null 
topological sector containing no universes at all to the sector with one universe of topology 
S3. The surface a = —00, \<f>\ < 00 can be thought of as a boundary between the two 
sectors. The probability flux is injected into superspace through this boundary (see Fig.l) 
and flows out of superspace through the remaining boundary (a — > — 00 with \<f>\ — > 00, 
or a — > +00). One could have thought that topology- changing transitions always occur 
through the boundaries of the corresponding superspace sectors. This was the point of view 
I adopted in my earlier formulation of the tunneling boundary condition 1 ^. I no longer 
believe this picture to be correct, but it may still be useful in cases when topology change 
is a semiclassical tunneling event. In this section I shall first review the motivation for the 
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old approach, then explain why I think it is not applicable in the general case, and finally 
discuss some alternative approaches to topology change. 

Tunneling amplitudes in quantum field theory are often evaluated semiclassically using 
the steepest descent approximation. One then finds that the path integral for the amplitude 
is dominated by a solution of Euclidean field equations, called the instanton. If topology 
change is a quantum tunneling event, one can similarly expect it to be represented by a 
smooth Euclidean manifold, Ai, interpolating between the initial configuration X4 and the 
final configuration, £2. The intermediate superspace configurations can be obtained as slices 
of M. and can be conveniently described using the concepts of Morse theory 47 . 

Consider a smooth real function f(x) on manifold M.. A point xq is called a critical 
point of / if d^f(xo) = 0. A critical point is called non-degenerate if det[d^d v f(xo)\ ^ 0. 
We shall call f(x) a Morse function if it has the following properties: (i) f(x) takes values 
between and 1, with f(x) = iff x G Si, and f(x) = 1 iff x G £2; (") a U critical 
points of / are in the interior of M. (that is, not on the boundary) and are non-degenerate. 
In a 2-d example of Fig. 2, the manifold M. is shown embedded in a 3-dimensional space, 
and the Morse function is given by the projection on the vertical axis. In this case, the 
saddle point P is a critical point of f(x). It can be shown that a Morse function can always 
be defined and that it always has some critical points if Si and S2 have different topology. 
We shall assume that f(x) is chosen so that it has the smallest possible number of critical 
points, that is, no more than dictated by topology. 

Slices of M. corresponding to superspace configurations can be obtained as surfaces of 
constant /. (Different choices of Morse function will, of course, give different slicings). 
These slices will have a smooth geometry, except the critical slices passing through critical 
points. With an appropriate choice of locally-Cartesian coordinates, the Morse function in 
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the vicinity of the critical point can be represented as 

d 

f(x) = J2^h (7-1) 
i=\ 

where d is the dimensionality of space (d — 3). The critical section, f(x) = 0, is a generalized 
cone. For d > 3 it has a curvature singularity of the form 48 



change is the "creation of universes from nothing", when the initial configuration is absent. 
A 2-dimensional illustration is shown in Fig. 3. Here, the critical slice is a single point. 
For near-critical slices, in d > 2 the curvature is again given by (7.2), where now r is the 
characteristic size of d-space. 

The idea of Ref.17 was that the boundary of superspace can be divided into regular and 
singular parts. The regular boundary includes only configurations which can be obtained 
as critical slices of smooth Euclidean manifolds. Such configurations correspond to transi- 
tions between different topological sectors. The remaining part of the boundary is called 
the singular boundary, and the outgoing-wave boundary condition is imposed only on that 
part. The boundary condition on regular boundary was supposed to enforce conservation of 
probability flux as it flows from one topological sector to another, but no specific form of 
the boundary condition was proposed. The overall picture was that the probability flux is 
injected into superspace through the boundary with the null sector, it then flows between dif- 
ferent topological sectors through the regular boundaries, and finally flows out of superspace 
through the singular boundary. 49 

As I mentioned earlier, I no longer think this picture can be valid in the general case. 
The main reason is that topology change does not necessarily occur between configurations 



R oc r 



-2 



(7.2) 
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at the boundaries of superspace sectors, but generally involves configurations in the interior 
of these sectors. It is true that, in order to change topology, one has to go through a 
singular 3-geometry. But, as we discussed in Section 6, superspace includes a very wide 
class of configurations, such as metrics with integrable curvature singularities and scalar 
fields with discontinuous derivatives. Note in particular that curvature singularities (7.2) on 
critical slices are integrable, and therefore the critical slices will generally lie in the interior 
of superspace. 

To give a specific example, consider creation of a wormhole in a universe having initially 
the topology of S3. The transition is then between the topological sectors S3 and Si x S2. 
The wormhole radius can be defined as r = (A TO j n /47r) 1 / 2 , where A m i n is the smallest cross- 
sectional area of the wormhole, and can be used as one of superspace variables. Since r has 
a semi-infinite range, r = is a superspace boundary in the sector Si x S2. On dimensional 
grounds, the curvature in the wormhole vicinity is ~ r , and the integral of 
does not diverge as r — > 0. The boundary at r = is therefore similar to what was called 
the regular boundary in Ref.17. On the other hand, configurations in the S3 sector "right 
before" topology change do not lie on any boundary. These configurations should only satisfy 
the continuity requirement: all matter fields should take the same values at the points that 
are about to be identified. 

An important example of topology change in lower dimensions is reconnection of inter- 
secting strings. At the classical level, this process plays a crucial role in the evolution of 
cosmic strings ^ (see Fig. 4). At the quantum level, it represents the elementary interaction 
vertex in fundamental string theories. A string loop can be thought of as a one-dimensional 
closed universe. The superspace configurations for the loop are given by the functions x^(a), 
where a is a parameter on the loop and the spacetime coordinates x^ play the role of world- 
sheet scalar fields. Topology change (loop splitting) can occur in configurations where the 
loop self-intersects, that is, when x^{a\) = x^{a2) for some o"i,o"2- These configurations are 
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not special in any other way and do not lie on superspace boundary. The configurations 
immediately after splitting have discontinuous derivatives of x^(a) at reconnection points. 
They are also legitimate superspace configurations and do not belong to a boundary. 51 

The conclusion is that topology-changing transitions affect not only superspace bound- 
ary, but can occur between points in the interior of different topological sectors. This has 
an important implication: in order to account for topology change, the Wheeler-DeWitt 
equation has to be modified. In the tradition of the subject, I would like to offer some 
speculations regarding the form of this modified equation. 

My suggestion is that the Wheeler-DeWitt operator Ti in (1.2) should be modified by 
adding an operator 5 that has matrix elements between different superspace sectors. The 
corresponding action can be written symbolically as 

[dh}i>*Hi) + / [dh] [dh'ty* (h)8(h, h')^(h') , (7.3) 



where the integration is taken over all superspace sectors and h stands for all superspace vari- 
ables. It seems reasonable to assume that topology change is a local process, then we should 
have 5{h\,h2) = unless hi and can be obtained from one another by changing topo- 
logical relations at a single point. The Wheeler-DeWitt equation for ipjqih) in topological 
sector iV is obtained by varying (7.3): 

HMh)+ f[dh']5 NN ,(h,h')1>N<ti) = 0. (7.4) 

The form of the operator S(h, h') is, of course, unknown. One can hope to gain some 
insight into it by studying lower-dimensional examples. In the case of strings, the topological 
sectors can be labelled by the number of disconnected loops, n, and 5 nn / has matrix elements 
with n = nil. An even simpler example is given by pointlike particles, which can be 
thought of as (0+l)-dimensional universes with spacetime coordinates playing the role of 
scalar fields and the Klein-Gordon operator V 2 + m 2 playing the role of the Wheeler-DeWitt 
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operator. Topology change corresponds to elementary particle interactions, like the one 
illustrated in Fig. 5 for a A0 3 theory. Here, two paticles merge into one, and the 5 operator 
should be proportional to ^-functions ensuring that the initial and final particles have the 
same coordinates at the moment of interaction. It would be interesting to develop the first- 
quantized formalism for particle interactions in the form (7.4) and verify its equivalence to 
a quantum field theory with non-linear interactions. The possibility of equivalence between 
a linear system of equations (7.4) and a non-linear field theory may seem rather unlikely. 
It is well known, however, that the full content of a perturbative quantum field theory can 
be expressed as an infinite set of linear relations between the Green's functions (Schwinger- 
Dyson equations). A similar representation has also been obtained in matrix models of 
two-dimensional quantum gravity 52 . 

To formulate the boundary conditions for the functions ip]^{h) in (7.4), we again divide 
the superspace boundary into singular and regular parts. The singular boundary includes 
configurations with \U\ — > oo and the null part of the boundary at a — > — oo (see Section 6). 
The functions ipjy should have only outgoing waves at the singular boundary. These waves 
carry the probability flux, 

Jn = K^N^^N ~ ^N^^n) » ( 7 - 5 ) 

out of superspace. The waves flowing into and out of the regular boundary correspond to 
transitions between topological sectors. In the example of the S\ x S2 — > S3 transition, the 
flux flowing into the regular part of the boundary at r = in the S\ x 5*2 sector reappears 
through the source term on the right-hand side of (7.4) in the S3 sector. The boundary 
condition at r = should enforce flux conservation between the two sectors. I will not 
attempt to write down a specific form of this boundary condition. 

Assuming that outgoing waves can be defined along the lines of the previous section 
and that the flux conservation condition is formulated, one can hope that the wave func- 
tion defined by eqs. (7.4) is equivalent to the one given by the path integral (1.3), where 
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the integration is performed over 4-manifolds of arbitrary topology. It is known that any 
Lorentzian metric interpolating between two compact spacelike surfaces of different topology 
must either be singular or contain closed timelike curves 53 . The singularities, however, can 
be very mild 54 , and there seems to be no reason for excluding the corresponding spacetimes 
from the path integral. If all metrics of finite action are included, this would be more than 
sufficient to permit Lorentzian topology change. 

8. Comments on Third Quantization 

It has often been argued !9, 20, 21, 22, 23, 24 ^at an adequate description of topology change 
can be given in the third- quantization approach, where the wave function ip is promoted 
to the status of a quantum field operator. Topology change is then accounted for by self- 
interaction of ip. For example, a ip^ interaction allows a parent universe, say of topology 
S3, to split into two daughter universes of the same topology. This is probably adequate for 
one-dimensional universes (strings), where topology is characterized simply by the occupa- 
tion number of closed loops. However, in higher dimensions the situation is not so simple. 
For two-dimensional universes, one would have to introduce an additional field creating and 
annihilating handles, while three-dimensional topologies have not yet been classified, and 
one may need to introduce an infinite number of fields and interaction types. It is not evi- 
dent, therefore, that third quantization offers any advantages in describing topology change, 
compared to the "first quantized" approaches like (7.4) or (1.4). 

I would also like to comment on the specfic implementation of the third quantization 
picture in simple minisuperspace models 19 ' 21 ' 24 . Without introducing non-linearity, the 
creation of universes in this approach is described in a manner similar to the description 
of particle creation in a time-varying external field. The idea is suggested by the fact that 
the Wheeler-DeWitt equation is similar to Klein-Gordon equation with the scale variable 
a playing the role of time and the superpotential U playing the role of a time-dependent 
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external potential. For a — > — oo the potential vanishes (see Sec. 6), and one can expand 
the filed operator ip into positive and negative-frequency modes, 

k 

with ipk( a ~ > — oo) oc exp(ia;^Q;) , a;^ > and creation and annihilation operators satisfying 
the usual commutation relations. The "in- vacuum" state, containing no universes at a — > 
— oo, would then be defined by a^|0 >{ n = 0, and single-universe states would be given by 
| A; >= a^|0 >. In the opposite limit of a — * +oo, one can similarly define a complete set of 
mode functions V'fejV'fcj such that ~~ * +oo) oc exp(iS) with dS/da > 0, and write 

^ = X^* + a *$P- (8-2) 
k 

The state containing no universes at a — > +oo is then |0 > ui with a^.|0 >out= 0, and 
single-universe states are |0 > ui- 

Since both sets of functions are complete, they must be linearly related to one another, 

4>k = ^(^kk'^k' + Pkkf^V) > ( 8 - 3 ) 
k' 

and eqs. (8.1), (8.2) then imply a linear relation between the creation and annihilaiton 
operators, 

fc' 

If the universal field ip is in the state |0 >j n containing no universes at a — > — oo, then the 
average number of universes in state k at a — > +co is generally non-zero and is given by 

< "ft >=< 0|5+S fc |0 > in = |/3 fcfc /| 2 . (8.5) 

fc' 

The suggestion in Refs. i9 > 2i ; 24 i s that < > should be interpreted as the number of 
universes created from nothing. I disagree with this interpretation for the reasons that I will 
now explain. 
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In the third quantization picture, there are no universes of vanishing size (a — > — oo), 
and as a grows, the number of universes increases and finally reaches its asymptotic value 
< h\. > at a — > +00. The universes are created at finite values of a, that is, with a finite size. 
This is drastically different from the creation-from-nothing picture, where the universes start 
at zero size and continuously evolve towards larger sizes, so that all the "creation" occurs at 
a — > —00. 

The origin of the difference between the two pictures is in the fact that the "time" a is 
not really a monotonic variable: the universes can both expand and contract. The positive- 
and negative-frequency mode functions -?/>fc and ^ correspond, respectively, to expanding 
and contracting universes. From this point of view, what is described in third quantization 
as creation of a pair of universes at some a = cxq, is simply a contracting universe that turns 
around and starts re-expanding at a = ckq. 

This can be illustrated using a (0+l)-dimensional example: pair creation in an external 
field. Following Feynman, antiparticles can be interpreted as particles travelling backwards 
in time, and pair creation corresponds to a particle trajectory like the one shown in Fig.6. 
The trajectory can be represented as rr^(r) with —00 < r < 00. Using the string theory 
language, r is a worldsheet time coordinate, and x^ are target space coordinates. For an 
observer riding on the particle, r is a suitable time coordinate and x^{r) is a set of interacting 
scalar fields. ^ The field x°(t) decreases with r at r — > —00 and grows at r — > +00. On the 
other hand, an external (e.g., human) observer, whose home is in the target space, will use 
x° as his time coordinate. 

In the third quantization picture, the variable a plays the role of target-space time, x°. It 
is not impossible that some super-human observer living in this target space will observe the 
creation of pairs of universes. 56 However, we are interested in what happens from the point 
of view of a worldsheet observer, living inside the universe and using the worldsheet time r. 
In any case, it appears that the process described by the third quantization formalism (8.1- 
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5) does not correspond to a topology- changing nucleation of the universe that the authors 
f 19,21,24 k a( j in m i n( j_ 

9. Conclusions 

The wave function of the universe ip can be obtained either by solving the Wheeler - De 
Witt equation with appropriate boundary conditions or by performing a path integral over 
an appropriate class of paths. Our discussion in this paper was focussed on the tunneling 
proposal for ip. Although little was proved, our discussion lead to several conjectures which 
will be briefly summarized here. 

In the path integral approach, the tunneling wave function is defined as a sum over 
Lorentzian 4-geometries interpolating between a vanishing 3-geometry (a point) and given 
3-geometry. The sum is, in general, performed over manifolds of arbitrary topology. I have 
argued that the wave function defined in this way should satisfy the outgoing - wave condition 
on a part of superspace boundary. 

Superspace can be divided into topological sectors, and part of its boundary can be 
thought of as the boundary between different sectors. We call it regular boundary. The 
rest of the boundary, which includes "incurably" singular configurations, is called singular 
boundary (see Sec. 6 for more details). The outgoing - wave condition should be satisfied 
only on the singular boundary. I have argued that the superpotential (6.4) of the Wheeler - 
De Witt equation either vanishes or diverges almost everywhere at this boundary and that 
this may enable one to give a precise definition of outgoing waves. 

If the topology of the universe is restricted to be that of a sphere, then the outgoing - wave 
boundary condition may be sufficient to determine the tunneling wave function. However, 
in the general case, this condition has to be supplemented by some boundary conditions at 
the regular boundary. 



33 

If topology change is allowed, then I have argued that it will occur not only through the 
boundaries between the superspace sectors, but will generally involve configurations in the 
superspace interior. This will result in a modification of the Wheeler - De Witt equation. A 
possible form of the modified equation is suggested in Sec. 7. 

Apart from the tunneling approach, I gave a critical discussion of the "generic" boundary 
condition (in Sec. 3) and of the third quantization picture (in Sec. 8). 

The tunneling approach to the wave function of the universe was motivated by the 
analogy with bubble nucleation and we may still gain important insights into the complicated 
issues of quantum cosmology by studying the wave function of the nucleating bubble. We 
may also learn a great deal from quantum gravity in two dimensions, which can be thought 
of as quantum cosmology of one-dimensional closed universes (strings), and even from the 
ordinary quantum field theory, in which the branching propagator lines in Feynman diagrams 
can be thought of as branching O-dimensional universes (particles). However, in pursuing 
these analogies, one should remember that in all these cases the observer is usually assumed 
to be in the target space, while in quantum cosmology the observer lives on the worldsheet. 
The relation between the wave functions of the universe (bubble, string, particle) obtained 
by these different observers is an intriguing problem for future research. 
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Figure Captions 

Fig.x The probability flow in the minisuperspace model (5.1). In this conformal diagram, 
a plays the role of time and the role of a spatial coordinate. 

Fig. 2 Topology change in two dimensions. The initial configuration £1 has topology 
S\ and the final configuration £2 has topology S\ © S\. The manifold M. interpolating 
between £1 and £2 is shown embedded in three-dimensional space. The Morse function 
f(x) is given by the projection on the vertical axis. The critical point P and the critical 
section f(x) = f(P) are indicated. 

Fig. 3 Creation of a two-dimensional universe from nothing. Here, the manifold M. has 
a single boundary E, and the critical section consists of a single point P. 

Fig. 4 A loop of string intersects itself and splits into two. Sharp angles (kinks) formed 
at the point of reconnection propagate around the "daughter" loops at the speed of light. 

Fig. 5 This c/> 3 interaction diagram corresponds to two scalar particles merging into one. 

Fig. 6 Feynman's picture of pair creation in external field. A particle travelling backwards 
in time from t = +00 turns around and travels back to t — +00. 
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